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ABSTRACT 

We investigate the anti-correlation between faint high redshift QSOs and low- 
redshift galaxy groups found by Boyle, Fong & Shanks (1988), on the assumption that 
it is caused by gravitational lensing of a flat QSO number count, rather than by dust 
in the galaxy groups, or any other systematic effect. Using an isothermal sphere lens 
model, the required velocity dispersion is a = 1286t^ km s -1 . With an isothermal 
sphere plus uniform density plane, the velocity dispersion is a = 11431^3 km s_1 ' 
while the plane density is S c = 0.081 ±0.032/i g cm" 2 . Both these values for the veloc- 
ity dispersion are considerably larger than the ~ 400 — 600 km s" 1 expected for poor 
clusters and groups and imply that the mass associated with such groups is ^4x larger 
than inferred from virial analyses. If due to lensing, this measurement clearly tends 
to favour high values of f2o- We demonstrate how an estimate of flo may be obtained, 
finding the relation Q = 1.3 (n/3 x 10~ 4 /i 3 Mpc" 3 )(r/l hr 1 Mpc)(er/1286 km s" 1 ) 2 
where r is the extent of the anti-correlation and n is the space density of groups. 
In the current data systematic errors in the determination of n and r may dominate 
this measurement, but this will be a potential route to estimating fig in improved 
galaxy-QSO datasets where these systematics can be better controlled. 

We have compared our result with that of Williams & Irwin (1998) who find a 
positive correlation between bright LBQS QSOs and APM galaxies. Because the QSO 
number counts are steeper at bright magnitudes, there is no contradiction between this 
result and our own. Indeed, adapting the lensing analysis of Williams and Irwin to 
our use of groups rather than galaxies, we find that there is good agreement between 
the amplitude of the positive cross-correlation found for the bright QSOs and the 
amplitude of the negative cross-correlation found for the faint QSOs. This analysis 
leads to a common estimate of SIqcts ^3 — 4. This, however, is significantly higher 
than indicated from several other analyses. Further tests of the accuracy of the galaxy- 
QSO cross-correlation results and thus their implications for f2o and as will soon be 
available from the new 2dF QSO catalogue. 
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1 INTRODUCTION 

Gravitational lensing by galaxies and clusters produces two 
different effects in QSO surveys. At bright magnitudes, 



where QSO counts are sleep, a positive correlation of QSOs 
and for eground galaxies or dttsters can be produced, as ob- 

jects intrinsically fainter than the magnitude li mit are ampli 



fled a nd hence artificially added to the sample ( Gott & Gunn 



1974). At fainter magnitudes, where the QSO number count 



slope is much flatter, it is the reduction of observed area 
behind the foreground lenses which dominates, producin g a 
deficit in the background QSO number count (Wu 1994). 

Here we interpret the faint QSO-galaxy group anti- 
correlation result of Boyle, Fong & Shanks (1988) (hereafter 
BFS88) (see also Shanks et al. 1983 and Boyle 1986) in terms 



of gravitational lensing. This result was first interpreted in 
terms of dust in foreground galaxies and clusters obscuring 
background QSOs. However, observations of galaxy groups 
and clust ers do not show significant amounts of dust (Fergu- 



son 1993) and the limits are at a level which make the dust 
hypothesis uncomfortable. Previously Rodrigues- Williams 
and Hogan (1994) have suggested that the anti-correlation 
result may be due to lensing and this is the avenue we shall 
pursue here. The results in this paper are based in part on 
those of Croom (1997). 

If correct, the lensing hypothesis would allow important 
constraints to be placed on cosmology and large-scale struc- 
ture. The deficit of QSOs near a group or cluster can be 
used to weigh that structure. This method has the advan- 
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tage over other mass estimates in that it allows a measure- 
ment of the absolute mass of the cluster, while other esti- 
mators such as the measurement of velocity dispersions and 
the observation of shear due to strong lensing ar e effectively 
measur ing t he gradient of the cluster potential (Broadhurst 
et al. ]|99e|). Other authors (e.g. Taylor et al. 1998) have 
looked for a deficit of galaxies behind foreground clusters 
to measure the lensing magnification. The advantage of us- 
ing QSOs over galaxies is that they are easier to distinguish 
as background objects and their redshift distribution is well 
known. Of course, there is the disadvantage that QSOs are 
rare objects and so cannot be used to examine the mass dis- 
tribution of individual clusters, however they can be used to 
investigate the properties of a distribution of clusters. 

Previous searches for QSO-galaxy correlations at 
brighter magnitudes have produced varying results with 
most showing the observational evidence for QSO-galaxy as- 
sociations, (see Table 1 of Wu 1994) but the statistical basis 
for most of the results was limited. Recently, Williams & Ir- 
win (1998) have found a strong positive c orrelation between 
~ 60 B < 18 LBQS flHewett et al. 1995[ ) z > 1 QSOs and 
APM galaxies. Below we shall compare their results with 
ours to check whether these two observations provide a con- 
sistent picture for the mass distribution in the Universe. 

Section 2 reviews the lensing model we use in this paper. 
In Section 3 we compare these models to the BFS88 data. In 
Section 4 we compare the results of Williams & Irwin with 
those of BFS88. We present our conclusions in Section 5. 



2 STATISTICAL GRAVITATIONAL LENSING 

We use two analytic mass profiles to fit the observed anti- 
correlation; the first, and simplest, of these being the single 
isothermal sphere (SIS), which gives a gravitational lensing 
amplification of 



A 



e > 0e, 



(i) 



(e.g. Wu 1994) where #e is the Einstein radius, the radius 
within which multiple images can occur. For the SIS case 
this is 



E = 4tt 



(§)'■ 



(2) 



where D s , D\ and D\ a are the angular diameter distances 
from the observer to the source, the observer to the lens 
and the lens to the source, respectively. In our second mass 
profile we add a uniform density plane to the isothermal 
profile (SIS+plane). This could be a good approximation to 
the effects of clustering and large scale structure (as pointed 
out by Wu et al., 1996), because a distribution of isother- 
mal spheres with an auto-correlation function of the form 



produces a uniform mass surface density. The 



globally measured auto-correlation function slope is ~ —1.8 
(Davis et al. 1988), which produces a sheet of matter which 



is uniform to better than 10% on the scales of interest. The 
amplification then becomes 



A 



1 



(3) 



„ - E /(1 - Sc/E crit ) (1 - E c /E crit ) 2 ' 

(e.g. Wu et al. 1996) where E c is the mass surface density 
in the plane and the critical surface density, E C rit, is 



E cr it — 



D s c 2 



D 1s DAtvG 



(4) 



Gravitational lensing can cause an over- or under- 
density of source objects near to the lens. The ratio of ob- 
served surface density to the true surface density (unlensed) 
is the enhancement factor, q, given by 



<? = 



N(<m + 2.5log(A)) 1 



N(< m) 



A 



(5) 



(Narayan 1989), where A is the amplification factor. N(< 



m) is the integrated number count of source objects brighter 
than magnitude m. We note that q depends on the source 
counts fainter than the limit of the survey. With a number 
count of the form N(< m) oc 10 am , we then find an angular 
cross-correlation function ujcq{&) that is described by 



wcq(0) = q- 1 = A 2 



1. 



(6) 



3 THE CORRELATION OF DURHAM/ AAT 
QSOS AND GALAXY GROUPS 

We l ook at the r esult from the Du rham /AAT UVX Sur- 
vey (Boyle 1986; Boyle et al. 1990) which shows an anti- 
correlation between UVX QSO candidates and galaxy 
groups (BFS88). This cross-correlation was carried out 
within 7 UKST fields, using COSMOS scans of photographic 



plate s. Spectroscopy of the UVX catalogue (Boyle et al 
1987) suggested that with a colour limit of u — b < —0.4 
there was ~ 55 per cent contamination by Galactic stars. 
In the BFS88 analysis the UVX criterion was tightened 
to u — b < —0.5, reducing contamination to 25 per cent 
while keeping 85 per cent of the QSOs. The UVX cata- 
logue was then split into two magnitude limited samples, 
17.9 < b < 19.9 and 17.9 < b < 20.65. The galaxy cata- 
logue consists of all galaxies to a limit of &j = 20.0 and the 
cluster sample was created using a 'friends-of-frien ds' algo- 
rithm (|Gott fc Turner 197?|; [Stevenson et al. 1988[). Groups 



of seven or more galaxies with density greater than 8 times 
the average for the field were classed as clusters, which 
amounted to 10 per cent of the total number of galaxies. 
BFS88 performed a cross-correlation analysis between the 
entire galaxy catalogue and the UVX sample but no signifi- 
cant correlation was found on any scale. Cross-correlation of 
cluster galaxies with the UVX catalogue resulted in negative 
correlations on scales < 10' for both samples, the brighter 
sample showing a marginally more negative clustering sig- 
nal. This can be interpreted as a decrease in contamination 
from smaller photometric errors in the brighter sample, but 
a second effect is that the QSO N(m) slope will be steeper 
at this brighter limit, thus a smaller anti-correlation might 
be expected. Given that our results are sensitive to the ex- 
act shape and position of the break, we restrict our analysis 
to the fainter sample, with the proviso that in new larger 
samples the anti-correlation as a function of magnitude will 
provide an important test of the gravitational lensing hy- 
pothesis. 



3.1 A comparison of lensing models and the data 

The effect of gravitational lensing is strongly dependent on 
the slope of the QSO number-magnitude relation at faint 
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1 1.5 
redshift 

Figure 1. Redshift distributions assumed for the galaxy and QSO 
populations in the lensing models described in the text. The nor- 
malizations of the two curves are arbitrary. 



magnitudes. We use the number counts from Boyle, Shanks 
& Peterson (1988) which give an asymptotic faint end slope 
of ~ 0.28. We have confirmed that this is a reasonable rep- 
resentation of the integral QSO number count at ~ 1 mag 
fainter than our magnitude limit, the region from where 
we expect amplified QSOs to come, by using the deeper 
data of Boyle, Jones & Shanks (1991). A flatter slope would 
clearly reduce the lensing mass required. The separation of 
observer, lens and source also affects the lensing amplifica- 
tion. To take this into account in our model, we integrate 
the known QSO redshift distribution over the effective range 
of the Durham/AAT survey (0.3 < z < 2.2). This gives us 
an effective lensing amplification for a particular lens mass. 
For the galaxies we assume the analytic form of N(z) given 
by Baugh & Efstathiou (1993): 



diV 
~dl 



oc z exp 



z (m) 



3/2 



(7) 



where z c = (0.016(fcj - 17.0) 15 +0.046)/1.412. This is shown 
integrated to 6j = 20 in Fig, jl] along w ith a polynomial fit to 
the QSO N(z) (|5hanks fc Boyle 1994[). The two populations 



occupy almost completely independent volumes, less than 
1% of the QSOs are at z < 0.3 while less than 0.5% of the 
galaxies are at z > 0.3. We assume an Oo = 1 cosmology 
throughout this analysis, but it should be noted that when 
the lensing mass is at low redshift (z ~ 0.1) cosmology has 
a relatively small effect as D s ~ Di s . 

We compare the SIS lens model (from Eqs. [I] and ^ 
to the cross-correlation result lucq{6) of the faint sample 
shown in Fig. |^. We have allowed for 25% contamination of 
the QSOs by randomly distributed stars. Using a minimum 
X 2 fit the velocity dispersion is a = 1286tgi km s _1 (reduced 
X 2 = 1.44). The dotted line in Fig. || shows this model. The 
SIS+plane model (Eqs. ^ and ^) is shown as the dashed line 
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Figure 2. UVX-cluster cross-correlation for the faint UVX sam- 
ple of BFS88 showing the best fit models; the dotted line shows 
the best fit SIS model, while the dashed line shows the SIS+plane 
model. The long-dashed line shows a lensing SIS model with 
a = 600 km s _1 . Also shown as the dot-dash line is the best 
fit 9~ 0,8 power law. 
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Figure 3. Confidence contours for the fit of the two component 
SIS+plane model to the cross-correlation result for the faint UVX 
sample, corrected for 25% stellar contamination in the QSO sur- 
vey. Contours are at A\ 2 = 1, 2, 3,4 (A% 2 = 1 and 4 correspond 
to 1 and 2 a errors respectively). 
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in Fig. ^| the best fit in this case has velocity dispersion 
a = 11431^53 km s _1 and the surface density in the plane 
of E c = 0.081 ± 0.032h g cm" 2 . The reduced x 2 for this fit 
is 1.07. The confidence levels for this fit are shown in Fig. || 

The values found are significantly larger than directly 
measured velocity dispersions of poor clusters and groups. 
However, because BFS88 correlate cluster members with 
UVX objects, they effectively weight each cluster by the 
number of member galaxies. Stevenson et al., (1988) find 
the fraction of clusters as a function of the number of mem- 
bers follows an approximate power-law with a slope of —2.2. 
From this we can calculate the mean cluster membership, n. 
Integrating the relationship between n = 7 and n = 50 gives 
n = 15. However, a member weighted mean gives n = 20. 
Thus, the BFS88 result is probing clusters which typically 
have ~ 20 members. The density on the sky of these clus- 
ters is ~ 0.8 deg -2 , which can be compared t o the den- 
sity of clusters in the APM Cluster Catalogue (Dalton et 
al. 199'jj) of ~ 0.2 deg -2 and the density of richness class 
or greater Abell clusters which is ~ 0.1 deg -2 (Abell et 
al. 198!fl). Thus an 'average' cluster used by BFS88 is sig- 
nificantly poorer than Abell richness clusters. The velocity 
dispersion that might be expected for clusters of this rich- 
ness is a ~ 400 - 600 km s" 1 (e.g. Ratcliffe et al. 1998). 
For comparison, a lensing model corresponding to a veloc- 
ity dispersion of 600 km s _1 is shown in Fig. the model 
is formally rejected at > 5er. It therefore appears that the 
masses implied for the galaxy groups from lensing are ~ 4 
times bigger than expected from virial analyses. 

Although the addition of the uniform plane to mimic 
the effects of clustering of clusters helps to improve the fit to 
<^QC Fig. ^ shows that this only reduces the velocity disper- 
sion of the clusters for high values of S c . Wu et al., (1996) 
find that the maximum mass likely to be associated with 
lenses from large-scale structure is ~ 0.01 — 0.02ft g cm" 2 , 
which assumes that the matter density in clusters is near the 
critical density (i.e. f2 c ius ~ 1). Values of E c in this range 
are only compatible with group velocity dispersions greater 
than 1000 km s" 1 (see Fig. |) 

We now demonstrate how the lensing estimates of av- 
erage group mass, via the velocity dispersion a, could be 
used to obtain a new estimate of Qq. Using the 0.8deg -2 
sky density of groups from above, we infer an approxi- 
mate space density of galaxy groups in the range n = 
2 — 4 x 10 _4 h 3 Mpc -3 . The lower value comes from inte- 
grating the proper volume to z=0.1 and assuming all groups 
are detected to this redshift, and on the basis of Fig. hi that 
this contains half the group sky density. The higher value 
comes from using the galaxy n(z) in Fig. [j] to derive the 
galaxy selection function. This is multiplied by the proper 
volume and the product is integrated to z=0.7, to give the 
effective volume from which the group space density is then 
derived. Multiplying by the estimated mass per group ob- 
tained by integrating the isothermal sphere profile out to a 
radius r leads to the estimate of Qq. We therefore find 



Q = 1.3 



0. 0003/i 3 Mpc- 



l/r 



iMpc) ( 



1286 km s" 1 



where r is now the extent of the anti-correlation, or the ef- 
fective extent of the isothermal sphere. We note that the 
dependence on ft cancels out and there is no dependence on 
any biasing parameter. The above scaling values represent 



our best estimate for n, r and a. The value for r is obtained 
from consideration of Fig. ^| where 8 ~ 10' corresponds to 
1 ft -1 Mpc. The errors on n and r are unfortunately likely to 
be dominated by systematic components, with each poten- 
tially varying by a factor of 2. This is due to the approximate 
methods used to determine both the surface density of clus- 
ters, the space density of clusters, and the limiting scale of 
the anti-correlation. We also note that the measurement of 
a is dependent on flo through the angular diameter distance 
terms in eq. ^, although at the redshifts considered this ef- 
fect is small (a ~ 25% difference in mass between flo = 1 
and fio = 0), and may currently be dominated by more 
serious systematic effects. More meaningful error estimates 
must await larger galaxy-QSO datasets, possibly with full 
redshift information for the galaxies as well as the QSOs, 
where the extent of the anti-correlation and the group den- 
sity can be better defined. 

A further final problem is that this analysis for Qo 
also assumes that the groups are physically real and not 
significantly contaminated by accidental line-of-sight over- 
densities on the sky, which Stevenson et al (1988) suggested 
was a possibility. Of course, even accidental over-densities 
will act as lenses and it is not yet clear how sensitive this 
estimate of fio is to this type of contamination. Again, in a 
survey which also has full galaxy redshift information, such 
as the forthcoming 2dF QSO/galaxy redshift survey, the ef- 
fects of spurious groups could be more easily checked. 



4 COMPARISON WITH THE LBQS-GALAXY 
CROSS-CORRELATION RESULT 

We now compare our conclusions with those of Williams 
& Irwin (1998) (henceforth WI98) who have found a strong 
positive correlation between APM galaxies and LBQS QSOs 
which is significant out to scales of ~ 60'. These authors 
find that their positive correlation is an order of magni- 
tude larger than that expected from a model with Qq — 0.3 
and galaxy bias of b ~ 1, based on a comparison of the 
galaxy-galaxy angular correlation function and the QSO- 
galaxy cross-correlation function. WI98 derive the relation: 



uj qg (8) ~ (2r/6)(2.5a - 1)wgg(0). 



,(8) 



(9) 



Here a is the slope of the QSO number counts, b is the 
galaxy bias, assumed to be constant as a function of scale 
and r is the optical depth of the lenses: 



(cdt/dz)(l + zf 



dz. 



(10) 



For a given value of fio, a bias value can therefore 
be found. This analysis can easily be applied to the faint 
QSO anti-correlation result of BFS88. We fit a power law 
with a slope of -0.8 to the auto-correlation of clusters 
measured by Stevenson et al. (1988), finding ucc(d) = 
(0.140 ± 0.053)6>" - 8 . We then also fit a -0.8 power law to 
the anti-correlation between QSOs and clusters, which we 
find to be w G q(0) = (-0.0071 ± 0.0059)6»-° 8 . With an as- 
sumed number count slope of 0.28 ± 0.02 these values then 
imply a value of rib = 0.085 ± 0.077. We assume Qo — 1 
and integrate Eq. llOj to z maK = 0.2, the redshift at which 
the N(z) relation described by Eq. m falls to half its peak 
value, this gives r = 0.021. We therefore find a bias value 
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of the clusters used in this analysis of be = 0.25 ± 0.23. If 
we use an Qo = 0.3 model, as used by WI98, then r and 
therefore be will fall by a factor of ~ 3. We should note here 
that the errors are large, and there is some uncertainty in 
this procedure; if we integrate r to where the N(z) drops to 
3/4 (z = 0.16) or 1/4 (z = 0.24) of its peak value we find 
r = 0.015 and 0.029 respectively. However even if we take 
the largest reasonable value of r, then 6c = 0.34 ± 0.031, 
which is still an order of magnitude lower than expected 
for clusters. Thus, in rough agreement with ba ~ 0.07 from 
WI98, we find a bias value estimated from statistical lensing 
which is an order of magnitude less than predicted by other 
methods. We also note that the WI98 result is consistent 
with the QSO-galaxy cross-correlation measured by BFS88, 
although BFS88 do not find a significant anti-correlation. 

Although our result appears to be consistent with WI98, 
they are both clearly significantly out of line with other cur- 
rent estimates of the combination of Oo and bias. The space 
density of galaxy clusters gives fig'Vs ^ 0.5 (Eke et al. 
1998) and dynamical estimates such as the measurement of 
redshift space distortions give similar values (e.g. Ratcliffc 
et al. 1998). We could possibly appeal to the scale depen- 
dence of bias to bring these different results into agreement, 
however, this would require an order of magnitude change 
in bias over a scale of ~ 10 h^ 1 Mpc. Taken at face value 
the above lensing result appears to suggest much more mass 
is present in the Universe than is detected from the distri- 
bution and motion of galaxies. 



5 DISCUSSION AND CONCLUSIONS 

BFS88 originally interpreted the UVX QSO-cluster anti- 
correlation as being due to absorption by dust present in 
clusters, the required amount of absorption being Ab ~ 0.2 
mag. Ferguson (1993) finds no evidence for any reddening 
due to dust in clusters, and the 90% upper limit on the red- 
dening is E{B — V) — 0.06. This upper limit is just consistent 
with the required absorption assuming Ab = 4.10E(B — V), 
and it is therefore still possible that lensing and absorption 
could both play a part in producing the anti-correlation re- 
sult. However, it is impossible for dust in groups to also pro- 
vide an explanation for the strong positive QSO-galaxy cor- 
relation found by Williams & Irwin and if their result is due 
to lensing then an anti-correlation is expected at faint QSO 
magnitudes which is comparable to that discussed here. If 
both results prove to be real, the simplest interpretation is 
that both are due to gravitational lensing. 

Assuming that the measured anti-correlation is due 
to gravitational lensing, we find that fitting an isothermal 
sphere model for the cluster potentials gives a larger than ex- 
pected velocity dispersion. Adding a uniform density plane 
to the mass profile does not significantly affect this con- 
clusion. These lensing mass estimates suggest cluster/group 
masses which are ~ 4 times larger than expected from virial 
theorem analyses. We discuss a potential method to deter- 
mine Qo from this type of mass estimate combined with a 
cluster/group space density measurement. We demonstrate 
this method with the current data, although an accurate 
measure of fio will have to wait for larger and better con- 
trolled galaxy-QSO dataset. 

We find consistency between the high Qo/b ~ 3 — 4 value 



implied by the strong positive QSO-galaxy cross-correlation 
seen at bright QSO magnitudes (Williams & Irwin 1998) 
and the negative QSO-galaxy cross-correlation seen at faint 
QSO magnitudes (BFS88), if lensing is assumed to cause 
both effects. Applying the method of Williams & Irwin to 
both these cross- correlation results gives Qoo~8 ~ 3—4 (where 
as ~ 1/&) and the inferred values of Qq''<js are therefore 6-8 
times higher than those inferred from arguments based on 
the space-density of rich clusters. 

Of course, it is still possible that some combination of 
systematic and random errors have contrived to produce the 
positive QSO-galaxy correlation seen in LBQS and the anti- 
correlation detected by BFS88. The importance of the above 
results suggests that it is vital to make further observational 
checks as to the reality of the QSO-galaxy cross-correlation 
signal. Fortunately, extended analyses of the above type 
will soon be possible with the completion of new large red- 
shift surveys such as the 2dF QSO Redshift Survey and 
the 2dF Galaxy Redshift Survey. These two samples with 
25000 QSOs and 250000 galaxies covering the same areas 
of sky should allow a definitive measurement of the cross- 
correlation function between background QSOs and galaxies 
at both bright and faint QSO magnitudes. 
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